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preface 


- 

For  most  applications  of  satellites  and  re-entering  space  vehicles, 
one  mst  be  concerned  with  the  accuracy  with  which  position  and  velocity 
can  be  determined  and  predicted  while  on  orbit  and  at  the  time  of  earth 
Igqpact.  Thlo  Mamorandum  desorlbee  a  confuting  program^  for  estimating. 

In  terms  of  confidence  regions,  the  on-orblt  and  Impact  errors  of  such 
vehicles, 

m  estimating  Impact  errors,  guidance  errors  are  olned  with 
orbital  prediction  errors.  The  analytically  determined  sensitivity 
coefficients  are  used  In  this  program  as  a  means  of  error  propagation. 
Their  es^resslon  as  functions  of  orbital  parameters  nay  make  them  useful 
for  other  purposes,  such  as  estimating  performance  reqpilrements  of 
treudclng  and  prediction  systems. 
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SUMMARY 

IFhls  Memorandum  describes  a  conipatlng  program  for  determining  errors 
In  position  and  velocity  on  a  satellite  orbit.  Error  coefficients  are 
commuted  from  anaJytlc  formulas.  These  may  be  used  In  the  further 
ccaq^tatlon  of  systematic  and  random  errors  In  the  prediction  of  satellite 
position  and  velocity.  The  computing  program  handles  the  propagation  of 
variance-covariance  and  the  determination  of  confidence  regions  for 
position  and  velocity  estimates. 
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a  semi-major  axis  of  elliptical  orbit 
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x,y,z;  z^,y^,z^,  etc. 


true  azuMnaly  (orbital  parameter)  at  time  t^ 
coordinates  associated  with  Input  variances 
coordinate  of  point  on  orbit  at  time  t 
coordinates  of  point  on  orbit  at  time  t^ 
velocity  components 


^  angle  between  horizontal  direction  and  velo¬ 
city  increment 

7  angle  between  horizontal  direction  and  total 
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,  A2>  variances  added  to  x,  y,  z  coagponents  of  velo- 

^  city  errors  due  to  guidance  system 


X  with  subscript  -  eigenvalue  of  covariance  matrix 

P  angle  between  orbltrary  reference  line  and  radius 
from  earth  center  to  vehicle  at  time  t 


0^  (same)  at  time  t^ 

statistical  parameter  associated  with  a  par¬ 
ticular  distribution  function 


n  angle  which  the  nodal  line  makes  with  the 
reference  direction,  generally  through  the 
point  of  Arles 

a  angle  between  the  nodal  line  and  the  radius 
vector  at  perigee 
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1.  BffTROlXJCTlOH 

In  first-order  error  propagation  in  any  system,  a  necessary  step 
is  the  determination  of  partial  derivatives  which  are  error  sensitivity 
coefficients.  When  it  is  feasible,  there  is  an  advantage  in  having 
these  coefficients  expressed  as  analytic  fxmctlons.  The  coefficients 
concerned  with  the  propagation  of  position  and  velocity  errors  into 
position  errors  (and  other  related  coefficients)  for  Keplerian  orbits 
were  given  in  an  earlier  RAID  paper.  For  the  present  Memorandm  this 
work  has  been  extended  to  Include  veloei'^  exTors  at  a  teniinal  xxxint 
as  well  as  position  errors  on  the  earth' s  surface  for  inqpact  tra¬ 
jectories* 

These  coefficients  are  used  in  the  propagation  of  variance*co« 
variance  for  position  and  velocity  errors  in  orbits.  The  resuxting 
variance- covariance  matrices  are  used  to  determine  confidence  regions 
for  position  and  velocity  errors  at  selected  points  on  an  orbit.  The 
introduction  of  guidance  errors  for  orbits  that  are  Impulsively  changed 
permits  an  assessment  of  errors  for  points  on  ajMwr  trajectory  in¬ 
cluding  errors  at  Impact  with  the  earth's  surface,  when  the  new  tra¬ 
jectory  Intersects  the  earth's  surface. 

ntis  Naaorandua  is  intended  to  provide  sufficient  Infoniation  for 
possible  future  users  of  the  computing  program  to  assemble  the  proper 
input  data  and  Interpret  the  output  data.  It  should  also  provide  the 
equations  and  background  information  for  re-progrannlDg  for  another 
conqniter. 
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II.  ERROR  SEHSrriVITIES 

The  error  seneitlvltieB  are  coefficients  In  the  first-order  error 
equations  and  are  oibtalned  hj  partial  differentiation  of  the  equations 
of  motion  expressed  in  a  particular  coordinate  system. 

The  secoDd-order  differential  equations  characterizing  tvo-hody 
notions  can  be  solred  to  give  the  position  and  velocity  of  each  body  as 
a  function  of  time.  Vfhen  these  are  solved  to  give  the  motion  of  an 
infinitesimally  smaller  body  moving  about  the  center  of  a  body  of  great 
mass,  the  usual  fonn  taken  is  that  of  parameters  describing  the  shape 
of  the  path  (a  conic  section)  and  an  equation  (generally  transcendental) 
relating  time  and  angular  position.  Vlhen  the  "total  energy"  is  negative, 
the  path  is  an  ellipse  end  time  and  angle  are  related  through  Kepler's 
equation. 

COOBDIlIMB  afgCBKB 

Figure  1  shows  polar  coordinates  in  the  plane  and  a  graphical  re¬ 
lationship  between  true  anomaly  v  and  eccentric  anomaly  E.  Figure  2 
gives  more  detailed  position  and  velocity  coordinates  in  the  planes  and 
Fig.  3  shows  the  three-dimensional  picture.  In  all  representations 
Pl(x^,yi,Zl)  is  the  point  where  observations  are  made,  and  p(z,y,s) 
is  the  point  for  whiOh  predictions  an  made.  In  path  prediction  from 
initial  position  and  velocity,  it  is  essential  to  note  that  perigee 
and  apogee  are  initially  undetermined,  requiring  angiilar  position,  P, 
to  be  measured  from  an  arbitrary  reference. 


Fig.  3 — Three  dimenslonol  picture  of  trajectory 
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OOTLIHE  of  DERIVATIOHS 

The  two  parameters  describing  the  elliptical  path  are  a,  the  semi- 
major  axis,  and  e,  the  eccentricity.  In  terms  of  Initial  position  and 
velocity  (see  Fig.  2)  these  can  be  expressed  as 

12  'l  * 

- ? - 

2  k 
1  *’l 


2  ""l  ""l 

e  -  1  -  -5  — 


■ysT 


With  zero  time  reference  at  perigee  or  perihelion ,  etc«,  the  re« 

V 

lation  between  time  and  angular  position  is 


7^ 


t  •  E  -  e  Bln 


kxiovnas  Kepler's  equation.  The  eccentric  anomaly  £  Is  related  to  true 
anomaly  v  by  cos  E  v  +  e 


,  an  ambiguous  eo^resslon  unless  ve  Incltide 


1.+  e  cosv 

o  <  E  <  *  for  o  <  V  <  * 

X  <  S  <  2i(  for  X  <  V  <  2* 


AuxlltSLTy  relations  used  are: 


r  .  ■■2  J 


i  +  a  cos  V 


k  e  slg  V 

Ul  -  e^) 


coe  E  > 


a  r 

ae 


a(l  -  e^)  >  r 
er 


cos  V 


6 


tan  7 


e  Bln  V 
1  +  e  COB  V 


COB  7  = 


1  +  e  COB  V 
2 

1  +  e  +  2e  COB  v 


For  earth  aatelllte  ozhlts  and  halUatlc  nlBBllea^  k  equals  K  equals 

c 

.07^36^7^  for  tine  in  odnTXtes,  and  distance  In  earth  radius  units  eq:uaJs 
3444  n  ai. 

Using  the  constraint 

k(t  -  t^)  =  a^/^  +  e(sln  Ej^  -  sin  E)j  ■  constant, 

partial  differentiation  yields  the  first-order  error  espresaiona  for 
errors  in  a  and  e,  and  for  position  and  velocity  in  the  plane  at  any  tlan 
t  or  for  any  true  anomaly  v. 

The  initial  position  and  velocity  vectors  detemlne  the  plane  defined 

• 

by  a  unit  vector  n^.  Errors  in  coniponoi;t8  of  this  vector  axe  detexadned* 
Front  these,  errors  in  position  and  velocity  perpendicular  to  the  plane  and 
errors  in  angles  i  and  R  defining  the  location  of  the  plane  follow. 


The  flrst-oorder  error  equations  for  which  the  coefficients  have  been 
determined  analytically  by  partial  differentiation  are  listed  below. 

of  these  partial  derivatives  Is  quite  arbitrary  and  grew  out  of 
work  which  extended  over  some  time. 
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dA 

a 

d*iH 

Kgdi,  H 

e 

-K^ 

^1 

+  Kj  dij^ 

de 

“H 

^1 

+  di^ 

.  dr^ 

dr 

-  K.  — 
'us  a 

^*S.5^1 

dr 

*  Kj 

^1 

^  *  *16  ^  ^  *  "^8 

■M*  -  ♦  *10  '»*l  »  *U  ^1  *  *12 

^*1  “  h2  ^1 

^  "  hs  ^*1  ^  *24  ^*1 

dl  -  Kj  +  K^  dyj^ 

^  -  ^S  ^*1  ^  *R 

dy  »  Kjj  dl  +  K^  dfl 
dy  -  Kp  d^j^  +  K^  dy^^ 

di  •  dTj^  +  Dg  dj6j^  +  D,  dr^ 

dv  ■  ^  df +  Hg  d|5j^  +  Hj  dTj^ 

d7  -  df ^  +  tfg  d^j^  +  dTj^ 

In  tb«  following  error  equatlong,  tbe  coefficients  evaluated  later 


are  for  tbe  particular  ease  of  earth  satellites  or  ballistic  stlssiles  ifben 
units  of  displacsBMnt  are  n  ■!>  units  of  velocity  are  ft/see. 
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dy  -  Kjj  +  Kp 

dz  «  Kq  dij^  +  Kg  dij^  +  Kj  dz^ 
d*  »  Ni  di^  +  »g  (1*1  +  ^1 

d^  -  Pj^  d^j^  +  Pg  dy^ 

dz  -  di^  +  Lg  dL^  +  L^  dz^^ 

dfi  -  dij^  +  Pg  dij^  +  F^  dZj^ 

dy  -  Qj^  dz^  +  Qg  dXj^  +  Qj  dZj^  (in  degrees) 

dR  ■  dij^  +  Tg  dl^  +  T^  dXj^  +  Tj^  dz^ 

dR  >  error  in  range  on  the  earth's  surface  (in  plane  of 
tarajectoxy). 

ERROR  COBPPICnBTS 

nte  ejqpressions  for  coefflciraits  in  the  first-order  error  eq^uations 
which  are  partial  derivatives  are  listed  below.  In  sons  cases,  these  co¬ 
efficients  are  given  as  functions  of  others  which  are  also  listed: 


r .  -Zl 

*2  “  k 


2a 

X 


a(l  -  e^) 


9 


-  e‘)  eln 


If  ya(l  ^  e  ) 

^  -  -2 


It  e‘ 


a(l  -  e2)  3 


e  +  e  COB  +  2  cos 


% 

H 


a(l  -  e^)  (e  +  cos  v^) 
- - - - 

e  (1  +  e  cos 

A 

a(l  -  e  )  (e  +  cos  v^) 
(l  +  a  COB 


^•hh^^hhk 

ho-hh6*Shl^^h9 

*  ^3  Kj^g  +  Kg 

r®  Bln  E  3  a^^  k(t  -  tj  a  sin  E 
Kj^3  ■  r  -  x*  sin  “  2  r 


hk 


COB  sin  E  ae(8ln  E  -  sin  E^^) 
- ^  —  + - 5 


r^  Bin  S 
hs  *  r  Bin  1!^ 


Bin  E 

'  ■  '■■  -  COB  B 


10 


1  +  e  cos  1  +  e  cos  v 
*S.6  “  e  sin  ”  e  sin  v 


^  r,  cos  V,  +  2ae 
r  cos  V  +  2ae  1 _ 1 _ 

*S.7  *  r  sin  V  "  sin 


-(1  +  o  cos  Vj^) 
sin 

1  •>■  e  cos  V 
*^9  “  er  sin  v 


*23  *  “ 

Kg^  -  -  8in(£0  +  Vj^) 


K  =  ^  V  K  =  .0096745519 

A  y 

>^0  "V 


Kb"' *12 


11 


■'j  ■  ''so  '^23 

*Sc  "  ^ 

Kl-*20''s 

»M  ‘  S 


K,.  -rKji, 

N)  "  '  *23 

Kp  .  -Kjo  r  8jii(*  -  T,) 

K-  —K^  r  Bln(»  -  *i)  -  *22  '  "  ’l' 


% 


*^3 

sin  1 


lU 


ynr 


2a 


h’^ih*  ^2^“  h 
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Dg  -  Kg  +  Cg  Kq 

+  e  coe  V 

^'^T7~ 


D=  - 


■in  V 


5  (1  +  +  2  e  cos  v) 


,  ®2 
h’T 


F  3 

3  *y 

1  +  e  co«  V 
■  •  sin  V 


,  2  ae  r  cos  V 
«2  “  r  sin  V 

1  •>•  e  COB  V 
®3  “  e  r  sin  V 

"3  -  Oi  *3  *  “s  ■%  *  ®3  S 

X 

Si 

S3  -  (B|,  K3  *  » 

■  *A  ’  7 

®2  "  S  “ 


.000161(281 

.000164231 

.016636$ 
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®3“*c 


Lj^  «  Bin  r  +  (1\  +  Dj  K^)  S  COB  r 

Dg  Bln  7  ■>•  (Dj^  Hg  ■!•  K^)  8  cob  7 


D,  Bln  7  +  (D.  H,  +  D_  K>)  i  cos  7 

- 


■1“ 


H, 


'^1“ 


Jg- 


'^2- 


1^  cot  7  -  (Dj^  +  Dj  K^)  I  Bin  7 

Dg  COB  7  -  (Dj^  Hg  +  D5  K^)  fi  Bln  7 

Dg  COB  7  -  (I\  Hg  +  Dj  K^)  4  Bin  7 

-  ^20^^  COB  A  V  +  r  Bln  A  v)  Av  -  v  -  v 

•  • 

K2g(rj6  sin  A  V  -  r  cos  A  v)  -  Kgj^(*^  cos  A  v  +  f  sin  A  v) 

Bln  7  +  (D4  +  D^  K^)  ft  cos  7 

Dg  Bin  7  +  (D^  Hg  +  D^  K^)  ft  cos  7 

Dg  Bin  7  +  (D^  Hg  +  D^  K^)  ft  cos  7 

^Hg.DjKJ 

Di^Hg.DjK^ 
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III.  PROPAGATION  OF  VARIANCE- CX)VARIANCE 

This  program  starts  with  a  variance- covariance  matrix  for  posi¬ 
tion  and  velocity  errors  at  some  point  on  a  nominal  orhlt.  The  orhlt  is 
specified  by  position  and  velocity  at  the  point  and/or  its  osculating 
Keplerian  parameters.  If  necessary  this  variance- covariance  matrix 
is  transformed  to  a  new  cooirdlnate  system  and  also  transfoaed  with 
respect  to  units.  Transformations  to  new  coordinate  systems  or  from 
one  column  vector  to  another  are  obtained  by  the  matrix  multiplication 

T 

C  «  A  B  A 

vhere^  for  exsn^le,  B  is  the  variance-covariance  matrix  associated  with 
column  vector  x,  and  C  is  the  variance-covariance  matrix  for  errors  in 
a  column  vector  y,  and  A  is  a  sensitivity  matrix 


Tbe  Input  variance- covariance  matrix  may  come  from  a  mraber  of 
sources,  such  as  a  dliKerential  coorreetion  routine  ehldl  has  been  used 
to  process  eustual  tracking  data  and  has  the  variance- covariance  matrix 
for  Initial  condition  error  estimates  as  a  hy-produet  of  the  oxhlt 
determination  process.  Ve  have  used  for  our  source  a  program  vhleh 
simulates  the  errors  In  radar  tracking  of  a  satellite  axid  computes 
the  statistical  parameters  associated  with  least  sq.uarsB  polynomial 
fitting  of  short  arcs  of  the  trajectory.  Ve  are  now  using  a  noie 
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(2) 

general  program'  ^  which  generataa  varlance<-covarianeo  matrices  result¬ 
ing  from  the  use  of  a  vide  variety  of  tracking  data  from  as  many  as 
12  different  trackers  in  orbit  determination. 

*niis  input  matrix  represents  the  variances  and  covariances  for 
errors  in  initial  conditions  consisting  of  three  components  of  posi¬ 
tion  and  three  of  velocity.  Since  it  is  convenient  to  use  a  coordi¬ 
nate  system  associated  with  the  plane  of  the  trajectory  for  the 
determination  of  error  sensitivities^  a  coordinate  transfonaation  of 
the  input  variance- covariance  matrix  is  usually  necessary.  Since  the 
sensitivity  coefficients  as  given  in  Section  II  are  functions  of  the 
ICeplerian  parameters  a,  e^  ii  and  v^^  it  is  also  necessary  to  com¬ 
pute  these  parameters  from  the  initial  conditions,  x  ,  y  ,  z  ,  x  , 

vv 


(XXMBDUCTB  TBMSFOEMnnaSl  TO  HARE  OF  TIUlJlCTGRr 

Given  the  initial  conditions  x,y,z,x,y,z  in  an  arbi- 

o^  o'  o'  •'o'  o 

trary  inertial  coordinate  system,  the  transfomiation  of  the  initial 
variance- covariance  matrix  B  to  the  coordinate  systan  associated 
with  the  plane  of  the  trajectory  as  given  in  Figs.  1,  2,  and  3  is  ob¬ 
tained  as 


C  = 


A  o 
o  A 


o 


where 


A 


“2  “3 

^  “2  “3 
Pi  Pg  P3 
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with  the  elements  of  A  determined  from 


“l  =  n2P3  -  n^pg 

“2  =•  Vl  ■  V3 
“3  =  V2  ■  “2^1 

THE  OSCUIATIBG  KEPLERIAM  PARAMETERS 

To  obtain  a,  e,  i,  m,  from  the  initial  conditions 

*0*  ^o*  K* 


17 


1 


k  =  .07^4*36574  (when  displacements 
are  in  earth  radius  units  and 
time  is  in  minutes) 


where 


To  determine  i  and  co>  cooipute 

y  z  -  y  z^ 
*^00  00 
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sin  1  =  +  ;  n^y 


0  <  i  <  l8o 


COS  i  =  n. 


Iz 


i  =  ten 


.1  * 


"iz 


sin  d 


*  -1^*  «o  v  * 


sin  i 
o 


ly  "  yp  °13C  °ly 


Sin  d  >  o  if  >  o 
sin  d  <  o  if  z^  <  o 


cos  d 


yp  °lx  -  °ly 

r  sin  i 
o 


Oft  a  d  . 

TRAHSFOBNAXIOII  OF  VARIANCE-COVARIAlfCE  MATRIX  FRCM  imriAL  FOIBT  TO 

aboMer  trajectory  POBff 

If  represents  the  variance-covariance  matrix  for  errors  In 
the  initial  point  and  0  the  sensitivity  matrix  of  partial  derivatives 
vhlch  relates  errors  for  another  orbital  time  to  the  initial  condition 
errors,  then  the  variance- covariance  matrix  for  position  and  velocity 
estimates  at  the  new  orbital  time  is  given  by 

H  -  0  0*^ 


i 

t 


19 

vhere  the  elements  of  0  are  determined  from  the  error  coefficient  ex¬ 
pressions  and  more  speclflccOly 


The  H  matrices  resulting  from  transformations  of  this  kind  repre¬ 
sent  the  erzor  situation  for  other  trajeetooy  points  in  terns  of 
variances  and  covariances.  For  certain  purposes  a  further  des- 
scription  in  terms  of  confidence  regions  is  deslrahle.  It  is  possible 
to  an  er  per  cent  confidence  region  (generally  an  ellipsoid)  for 

position  errors,  or  another  for  velocity  errors.  This  Is  the  region  In 
which  estlaates  would  fall  er  per  cent  of  the  time  if  the  experiment  wore 
repeated  a  very  large  nosiber  of  tlxMS.  If  we  consider  a  partitioning 
of  the  E  matrix  Into 


then  represents  the  variance-covariance  auitrlx  for  position  errors 
and  the  varlanee-ecvarlance  matrix  for  velocity  errors. 
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The  quadratic  form  deflnlag  a  cmfldence  allipsold  for  position 

(7) 

errors  is  given  by'"" 


T 

X*  J, 


z 


2  2 

where  is  the  1-a  level  of  the  X  distribution  for  three  degrees 

o 

of  freedom.  For  or  «  ^  per  cent)  x  **  2.366  and  for  or  ■  9?  per  cent, 
2 

X  ■  7.619.  Accordingly,  the  eigenvalues  of  the  matrix  determine 
the  slse  of  the  seal-axes  of  the  confidence  ellipsoid  and  the  eigen¬ 
vectors  and/or  the  associated  rotation  matrix  determines  the  relative 
orientation  of  the  confidence  eUlpsold.  If  the  elganvalnes  are  re¬ 
spectively  X^,  X2*  X^  then  the  semi-major  axes  are  given  by 


dg  =72.366  Xg  B  =  1,2,3 

for  the  90  per  cent  confidence  ellipsoid,  and 

dg  »  Xg  a  m  1,2,3 

for  the  91$  per  cent  confidence  ellipsoid. 

Operating  in  an  identical  manner  with  determines  the  con¬ 
fidence  region  for  velocity  errors.  The  slx-dlaenslonal  confidence 
region  (hyperelllpsold }  for  the  combined  position  and  velocity  estimates 
is  obtained  in  an  analogous  manner,  with 

d„  =  75.3W  8  -  1,..  6 

O  O 

for  the  90  per  cent  confidence  region, 

and  j - 

d„  =  J12.992  X„  s  .  1,  ..  6 

B  »  O 

for  the  9?  per  cent  confidence  region. 
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aCORPORATIOW  OF  GUIDANCE  ERRORS 

When  a  trajectory  is  changed  by  Impulsive  velocity  components, 
errors  due  to  the  guidance  system  may  be  Introduced*  If  there  Is  no 
correlation  vlth  prediction  errors,  the  variances  in  toqiulslve  velocity 
components  are  simply  added  to  those  due  to  the  prediction  process* 

In  gen  ral  vhen  correlation  exists,  the  guidance  errors  ere  Incorporated 
by  transforming  to  covariance  matrix  for  a  nine  element  vector,  in« 
troduclng  guidance  error  variances,  and  then  transforming  back  to  a 
6x6  variance- covariance  matrix.  The  transfoxmatlon  which  we  use 
recognizes  the  possibility  of  a  relationship  between  the  predicted 
position  error  In  the  plane  and  in-plaae  velocity  component  error  when 
a  stellar  referenced  stabilized  platform  Is  used. 

The  transformation  required  Is: 

T 

M  «  L  H  L 

where  H  Is  the  variance-covariance  matrix  for  errors  at  the  trajectory 
point  before  the  impulsive  velocity  Increment  Is  added  and 
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-  V.  Bin  P 

-  _ _ -E - 

^71  r 

-  V  cos  p 

-<91 - V- 

See  Fig.  U  for  definition  at  P  and  V..  Iltals  Is  for  the  stellar  refer- 

p 

enced  platfom  and  applies  to  circular  orULta.  For  Ideal  eoapansatlon, 

p  and  V.  have  particular  values.  For  other  eases,  other  knosn  oorxela- 

P 

tions  could  be  Introduced  In  an  analogous  wanner. 

It  has  been  shown  by  Frick'  '  and  others  that  the  relationship 
between  the  impulsive  velocity  Increment  and  the  In-place  predicted 
position  error  for  circular  orbits  can  result  in  ccopensatlon  of  In 
plane  position  error  to  first  order  for  a  particular  velocity  in¬ 
crement  when  the  range  to  impact  is  fixed.  Figure  4  shows  sche¬ 
matically  how  this  is  accomplished.  This  range  error  compensation 
do 

occurs  vrtien  ^nd  i;he  particular  values  of  p  and  are  given 

by  solving  first  for  the  required  velocity  cooponents  as  follows, 


r^^tr^^  -  cos  AV)  x 


1  -  COB  dV  ;  2 

+  — 5 -  S 

r  ^ 


.  (1  -  cos  AV?^ -  ^  Q  ^ 

(r^  -  cos  AV) 


pnd 


1  (1  -  COB  AV)  1 

r,  sin  aV  a 


r^  -  cos  AV  I  . 
"sin  aV - )  ®x 
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vhsn  r  Is  the  ndiim  of  the  Initial  ocbit  saA  AY  is  the  nhen#t  in  tne 

anomaly  Identlcad  vlth  the  angiilar  range  to  liqpact.  Another  interestixig 

do 

case  is  the  minimum  ijopulBe  path  occurlng  vhen  ^  *  0*  Ihe  value  of 


and  ^  can  be  obtained  by  solving 
P 


To  Incoi^orate  guidance  errors  Into  the  M  matrix  we  add  variance 


terms  to  the  diagonal  elements  in  the  lover  right  hand  comer,  thus 

becomes 

becomes 


and 


^9 


becoMs 


^9^ 


If  coveurlances  in  the  guidance  errors  are  appreciable  and  known,  these 
may  be  added  alao.  The  rasulting  matrix  is  transferred  to 
a  6  X  6  by 


P  -  H 


with 


0 

1 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 


0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

1 

0 


0 

0 

0 

0 

0 

1 


0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

1 

0 


nopMuaiXM  or  yaeuahcb  n  aw  guuicioRr 


Since  the  sensitivity  coefficients  are  eiqpreaaed  anslyt^sUy  in 
terms  of  the  osculating  o  oital  parameters,  these  parameters  for  the 
new  orbit  must  be  determined.  These  follow  from 
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V  ■  +  Av 
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The  value  of  v  does  not  change  and  u)  is  arbitrary  for  the  pur¬ 
pose.  These  new  trajectory  parameters  now  constitute  the  input  to 
a  new  computation  of  error  sensitivity  coefficients.  If  one  is  con¬ 
cerned  about  the  variance- covariance,  for  errors  or  confidence  regions 
for  a  general  point  on  the  new  trajectory  the  procedure  is  identical 
with  that  described  cn  pages  19-20.  Hc^rever,  If  the  errors  In  a  tan¬ 
gent  plane  at  the  earth's  stirface  at  the  point  of  intersection  or 
impact  are  to  be  described,  the  transformation  requires  a  different 
set  of  sensitivity  coefficients.  The  2x2  varisnce»covarlance  matrix 
for  errors  at  the  impact  point  is  given  by 


R  =  Q  P  Q 

where  Q  is  the  sensitivity  matrix 

T, 


Q  s 


T-  0  T. 
3  A 


Ky  0 


and  the  sensitivity  coefficients  are  appropriate  for  the  new  trajectory 
and  the  impact  ix>int.  A  confidence  region  which  is  now  an  ellipse  is 
defined  by  the  quadratic  form 


R 


X 

y 


with 


=  1.386  for  a  •  50  per  cent  confidence 
X^^^  s  5.99  for  a  «  95  per  cent  confidence 


The  eigenvalues  and/or  eigenvectors  for  the  R  matrix  then  determine 
the  size  azui  direction  of  the  a€Bl-»axee  of  the  confidence  elUpee. 
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IV.  LAYOOI  CP  THE  FROQRAN 


The  routine  vms  coded  using  FAP  and  FORTRAN  for  the  IBM  7090 

cofliputer.  It  contains  the  following: 

Hand-coded  subroutines; 

PAST7  -  computes  the  orbit  change 

KEP  -  computes  error  coefficients 

A£I  -  conqputes  a,  e,  i^ 

NATNPYI  coaQ>ute8  product  of  two  matrloes  (6x6) 
NAINPBj  *  and  (9x9)  >  respectively 

STEP2  -  pre (post) -multiplier  of  the  Input 
covariance  matrix 

SCALE  -  scales  a  matrix  to  avoid  overflow  In  EIOSN 

ARCSiA  _  conq^utes  COS"^  €U3d  SIN"^  from  Hastings 

ARCCO^  ~  approximation 

AVQ6  -  averages  elements  of  a  real  syssetrlc  matrix 

XERA  -  checks  for  BEOIN  flag 

RTSXD  -  cooqputes  ^  and  Vg  on  ox^lt  change  (see 

options  for  thls^pcurt) 


Library  routines: 


SHARE  EIGEN  -  confutes  elgenvalxies  and  eigenvectors 
RAND  X006  -  tan*^  of  double  argument 

sortpI 

COSP  \  -  standard  FORTRAN  library 

SINP  J 

plus  the  master  routine 
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Step  1: 


Step  2: 


Step  5: 


Step  4: 


PROGRAM  DESCRIPTION 


Start  with  a  given  covariance  matrix^  B(6x6),  that 
gives  the  errors  in  the  Initial  conditions  of  a 
nominal  orbit  specified  by  ,  y^  ,  ,  x^  ,  y^  , 

These  conditions  are  combined  using  a  subroutine  called 
AEI  to  give  the  alternate  specification  of  the  orbit 
by  a,  e  ,  i  ,  a>  ,  v^^  .  The  eigenvalues  of  the 

matrix  B  are  found  and  printed*  If  these  not 
all  positive^  the  matrix  is  not  meaningful  and  the 
program  will  halt  later  on* 

A  transformation  matrix  A^  is  found  (subroutine  STEP2) 
such  that  the  given  covariance  matrix  is  tremsformed 
into  a  coordinate  system  associated  with  the  plane  of 
the  trajectory* 

This  new  matrix  is  now  called  C  where 

T 

C  «  A^BA^  . 

Convert  the  units  of  C  to  nautical  miles  €uid  ft/sec. 
to  give  the  matrix  C^^;  i.e*^ 


if 


then 


where  k. 


D  E 
E  F 

Dxk.  Sxk^ 

X  c 

Exk^  Fxk. 


3444  X  20.926  .  , 


20.926  X  10*^ 
- - 


N« 
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Step  Input  the  hv's  (diange  In  true  anonaly)  to  be  considered* 
Compute  the  error  coefficients  for  the  orbit  using 
subroutine  KHP.  Ihe  formulas  are  contedned  in 
Sec.  II  •  Although  all  of  these  coefficients  are  not 
used  in  the  program,  they  may  be  printed  out  at  the 
option  of  the  user  by  setting  the  correct  value 
of  KPRINT  (see  input  requirements).  The  sensitivity 
matrix  G  is  computed  using  the  proper  error 

coefficients.  Ihe  matrix  G  gives  the  proipa^ted  errors  In 

X,  y,  z,  X,  y,  z  ,  for  a  given  v  . 

Step  6:  For  each  G  compute 

T 

H  •  GC^G 

for  the  errors  in  position  and  velocity. 

Step  7:  Various  submatrices  of  H  are  then  used  to  give 

confidence  ellipsoids  for  which  semi-axes,  angles  of 
rotation,  eigenvalues  and  eigenvectors  are  confuted 
aind  printed. 

Itois  completes  the  first  part  of  the  program.  If  one 
wishes  to  go  on  to  a  new  trajectory,  the  value  of  KFLAG  is 
appropriately  set  on  input  and  computation  proceeds. 

Step  8:  There  are  three  options  for  input  of  data  on  the 


new  trajectory. 
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Step  9: 


Step  10; 


a.  specify  p  ,  ,  Av  ,  KLM  «  1  . 

b.  mlDlmum  inqpulse  path^ 

KLM  ■  2,  specify  A^  ,  A^  ,  A^  ,  Av  ,  machine 

confutes  ^  and  V^« 

P 

c.  range  conqpensatlon  path, 

KLM  »  3,  specify  A^  ,  ,  A^  ,  Av  ,  machine 

conqputes  ^  and  V.  . 

p 

Then,  having  a  value  of  ^  and  ,  an  L(6x9)  matrix 

P 

is  determined  so  that  the  M(9x9)  matrix  can  be  found 
M  -  LHL^  . 

is  formed  by  changing  three  elements  of  the  N  matrix 
as  follows: 


Vr  *°  *77**1 
1^  to  ^  ♦  ig 

■99  “  *99**5 

VtMZO  0^  ,  ^2  '  iiq>uts.  (Not  eumulAtlTc) 

A  new  transformation,  matrix  N(6x9)  ,  la  found  so  that 
the  following  transformation  can  be  made: 


P  - 

idiere 

1  0  0  0  0 

0  10  0  0 

,.00100 

0  0  0  1  0 

0  0  0  0  1 

0  0  0  0  0 


0 

0 

0 

0 

0 

1 


0  0  0 
0  0  0 
0  0  0 

10  0 
0  10 
0  0  1 
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Step  11: 

Step  12: 


Step  13: 

KPRINT  - 

KFLAG  < 
> 

KLM  -  1 

-  2 

-  3 

KMI 


New  values  of  a  ,  e  ,  1  ,  vtt  ,  v^  ,  and  new  values 

of  the  error  coefficients  (which  will  he  printed  If 
KPRINT  IS  SET  ^  O)  are  found. 

A  new  sensitivity  isatrlx  Q(6x2)  Is  coisputed  giving  the 
errors  In  range  and  y,  and  a  covariance  aatzlx  R(2x2) 

Is  determined  for  these  errors, 

R  -  QPq"^  . 

The  confidence  ellipsoid  seml-eoces  are  coiqputed  and  printed, 
vhlch  coaqpletes  the  problem* 

OPTIONS 

0  No  error  coefficient  printout* 

0  Error  coefficients  printed. 

0  Change  to  new  tre^)ectory  (ies)« 

0  Go  to  nev  case  -  no  change  of  trajectory. 

Specify  p  ,  Vp  ,  ,  Ag  ,  ,  on  new  trajectory. 

Specify  A^  ,  A^  ,  A^  ,  A^  .  coaqpute  ^  , 

minimum  isqpulse  path,  p  <  0* 

Specify  ,  Ag  ,  Aj  ,  A^  ;  coi^ute  P  ,  ,  range 

compensation  path.  P  >  0* 

0  Condensation  desired.  Omly  raaecsable  for  P  >  0  tnJeetoKlee* 
0  No  condensation  desired. 
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1- 


INPUT 


Card  No. 

Description 

Fomiat 

1 

U  0  Error  eoaffieisnts 

I  2 

KPRINT/  printed 

0  Error  coefficients 

not  printed 

2 

x^  ,  y^  ,  ,  (Earth  radius  units) 

SElB.d 

3 

\  radius  vnits/mln) 

3E16.8 

k^l^ 

Covariance  matrix,  B;  diagonals  are 

3E18.8 

varluteea  (or  standard  errors) |  off 
diagonals  are  cross  product  terns;  l.s., 


16 


Products 
of  the. 
units 


Products  of 
these  units 


.JCfflP/Bln.) 


-  Count  of  the  number  of  dr's 
(true  snoBMlles)  to  be  input. 


17 

e 

e 

17+l^UltT5- 


1 


17+1^^U1IT5 


dr's  in 


KFLAO 


1  par  card 


Change  trajectory. 

Don't  change  trajectory;  go 
to  new  case.  Next  card  is  the 


Ii»^X;0UNT  card. 
Ut  L  -  17+K0U1PP5H 


112 


F12.0 

112 


L  KM 


2112 


K^UIfF  *  Huiibtr  of  nw  trajoctorlos 
to  be  coaputed. 

■  0  cosgtensatlon  desired 
/  0  no  eoapensation  desired 


1+1 


KIM  -  As  explained  in  vrltev^  >  step  6 


I  12 


Card  No, 
L  +  2 


L  +  3 


LfigdunT 
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Description  Fonnat 

^1  >  Ag  /  Aj  /  A^  6EX2«8 

d-  ^  9  A-  ^  Elements  to  add  to 

M  aatrlx  to  get  aatzlx. 

-  Quantity  to  add  to  to  get  V. 

if  KLM  -  1;  P,  V-  6B12.8 

P 

if  KIM  =■  2  or  3  next  card  like  L  +  1 

Repeat  cards  L  1,  L  +  2  (and  L  3  if  KIM  >1) 
l^UR-l  tiaes. 

"BEQIH"  punched  in  cols.  1>^  •  Signifies  the  a6 
start  of  a  nev  case. 
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OUTPUT 


FIRST  PJCE 

1.  Characteristic  roots  of  B,  the  input  covariance  matrix. 
These  must  all  be  positive  or  B  is  not  a  valid  matrix 
for  this  problem.  The  first  three  are  in  ERU^  and  the 
last  three  in  (ERU/mln)^  . 

2.  The  initial  conditions  x,  y,  z,  x,  y,  z  in  ERU  and 
ERU/min,  respectively. 

3.  The  Keplerlan  paursMters  of  the  orbit. 

a  in  ERU 

e  non-dimensional 


k. 


5* 


The  input  coveurlance  matrix. 

p 

The  upper  left** band  comer  (3x3)  In  (XRtJ) 

p 

The  lover  rl^t-hand  comer  (5x3)  in  (ERU/mln) 

The  v^pper  right  and  lover  left-hand  comers  are 
combinations  of  these  units. 

Transformation  matrix,  . 

Non-dimensional,  takes  B  from  the  initial  reference 
system  to  a  coordinate  system  associated  vlth  the 
oxhital  plane. 
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SECOND  PAGE 

T 

!•  Transformed  matrix  C,  ■  AlBAl 
Units  are  the  same  as  B. 

2.  Matrix. 

p  P 

Units  changed  from  (EHU)  to  (nautical  mllea)  and 
(ERU/min)^  to  (Ft/ sec) ^  . 

3.  True  anomaly^  v,  change  in  true  anomaly^ 

The  angle  the  v^lcle  moves  from  perigee. 

U.  0  Matrix;  error  sensitivity  matrix, 
a.  First  row;  error  In  x  . 
h.  Second  row;  error  In  y  . 

c.  Third  row;  error  In  z  . 

d.  Fourth  row;  error  In  x  . 

e.  Fifth  row;  error  In  y  . 

r.  Sixth  row;  error  In  z  • 

T 

H  Matrix  ■  GCG  ;  transformed  matrix  gives  error  in 
position  and  velocity  at  a  new  point  on  the  trajectory 
specified  by  Av«  Same  units  as  • 

IKyiPE:  If  It  Is  desired  to  print  out  all  the  error 

coefficients,  they  will  be  printed  after  Step  2 
above,  and  Step  3  will  begin  a  new  page.  If  there 
Is  more  than  one  v  ,  each  new  v  will  begin  a 


new  page. 
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_ /Number  of  Av’s  +  2  If  error  coefficients  printed. 
- =S-  *  ^umber  of  Ar's  +  1  If  no  error  coefficients  printed. 

1.  (n)  (n  ■  1,  ...  number  of  Av's) 

Upper  left-hand  comer  of  H  . 

2 

Error  sltxiatlon  In  position  (n  mi)  e 

2.  Hoots  of  . 

3f  Seml*axes  of  confidence  ellipsoid^ 

l.e.,  where  you  would  expect  to  find  the  object 

!K>  per  cent  end  99  par  cant  of  the  tlaa^  raqpactival^e 

U.  L-Natrix  -  Eigenvectors  of  • 

Rotation  matrix  to  give  new  coordinate  system 
with  no  correlation. 

a  ,  ^  ^  7  -  angles  relating  new  region  to  the  old. 


N»2nd  Page 

Same  as  above  but  for  p  the  lover  right-hand 

comer  of  H. 

N»5rd  Page 

Same  as  above  but  for  total  H  matrix,  Nov  you  have  a 
confidence  hyperellipsoid.  Step  ^  is  not  done. 

NOTE;  These  three  pages  are  repeated  n  times  for  the  n 
values  of  v  . 


N»3n»lst  Page, 

!•  Vp  ,  P  ,  ,  Ag  ,  Aj  ,  Av  . 

P  «  angle  (in  degrees)  between  the  velocity  increment 
and  the  original  velocity  vector. 

Vp  •  velocity  increment  (ft/sec) 
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2. 


3. 


4. 


5. 


Ap>*  Erix>r  variances  for  guidance  (ft/sac) 


Av  *  Change  in  true  anomaly  on  nev  trajectory 
(in  degrees). 

Parameters  associated  vitb  nev  and  old  orbits  as 
labeled: 

a  ;  e  -  same  units  as  Step  5,  page  1. 

,  u)  ,  i  -  radians. 

Q(k^k)  matrix  -  another  sensitivity  matrix. 

First  row;  errors  in  r  . 

Second  row;  errors  in  y  . 

R(k,k)  -  covariance  matrix  for  r  and  y  errors. 

Semi-eoces  for  50  par  cant  and  95  par  cant  eonfldanca  as  bafcra* 


NOTE:  If  K  represents  the  number  of  nev  trajectories  to 

be  considered,  the  above  five  st^s  vill  be  output  K 
times.  As  before,  if  all  the  error  coefficients  are 
to  be  printed,  they  vill  be  printed  after  Step  2  and 
Step  3  vill  begin  a  nev  page* 
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NOTE-  ONE  MAY  OBTAIN  THE  OPTION  (DOHED  LINES)  BY 
REMOVING  THE  C  FROM  THE  APPROPRIATE  COMMENT 
CARDS  AND  INSERTING  C’s  ON  OTHERS 


GUESS  COMPUTE 
aiNCTION  + 
DEKtVATIVE 
FOKM  2 


Sk^  «  5x 
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SUBROUTINE  PAST  7 
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VI.  PRESEMT  LIMIIATIONS  AMD  POSSIBLE  MDDIfICATlOHS  OF  THE  PROCfflAM 

The  sensitivity  coefficients  are  here  determined  for  a  theoret¬ 
ical  Kepler ian  orbit.  Even  for  cases  of  moderately  high  dreg  the 
error  in  these  coefficients  due  to  this  assumption  will  not  be  serious. 
Hovever,  the  error  in  predicted  position  and  velocity  components  due 
to  error  in  estimation  of  drag  force  represents  an  additional  error 
component  that  is  not  Included  in  this  program.  An  effect  vhlch  la 
external  to  this  program  is  the  validity  of  the  input  covariance 
matrix.  This  program  does  check  on  the  requirement  of  positive  de¬ 
finiteness  of  the  input  covariance  matrix  for  physical  realizability, 
but  otherwise  imposes  no  restrictions. 

In  programning  prediction  intervals,  exact  multiples  of  180°  must 
be  avoided  since  the  sine  of  v  and  v^  occurs  in  the  denoeiinator  of 
certain  expressions,  rendering  them  Indetezninate.  The  error  In  longi¬ 
tude  of  the  node  also  becomes  Infinite  for  zero  inclination  angle. 

The  case  for  zero  eccentricity  or  exactly  circular  orbits  is  also 
avoided  in  the  program  by  eccentricities  never  less  than  .001 

in  computing  error  coefficients.  In  practice  this  gives  values  for 
error  coefficients  sufficiently  close  to  those  for  a  clreulsr  orbit. 

It  is  obvious  that  a  more  sophisticated  statement  of  guidance 
errors  including  their  covariances  could  be  incorporated  into  the 
present  program  without  much  difficulty.  A  major  modification  of  the 
program  would  be  required  to  cosqpute  error  sensitivity  coefficients  for 
the  high  drag  re-entry  case  since  this  would  involve  integration  of  the 
equations  of  motion  and  a  numerical  determination  of  partial  derivatives. 
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